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Abstract: We consider the growth of cosmological perturbations in modified gravity 
models where a scalar field mediates a non-universal Yukawa force between different matter 
species. The growth of the density contrast is altered for scales below the Compton wave- 
length of the scalar field. As the universe expands, the Compton wave-length varies in time 
in such a way that scales which were outside the range of the scalar field force may feel it 
at a lower redshift. In this case, both the exponent 7 measuring the growth of Cold Dark 
Matter perturbations and the slip function representing the ratio of the two Newtonian 
potentials ijj and <j> may differ from their values in General Relativity at low redshift. 
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1. Introduction 

Gravity has been precisely tested in the solar system with tight constraints on fifth forces 
and violations of the equivalence principle. Its non-linear regime has even been probed 
in the binary pulsar systems (for a review, see e.g. [1,2]). So far no significant deviation 
from general relativity has been detected. On the other hand, the discovery of cosmic 
acceleration in 1998 [3, 4] has prompted a rich and renewed bout of activity on models 
where gravity could be modified (for an overview, see e.g. [5] and references therein). This 
has been triggered by the difficulty of explaining cosmic acceleration with dark energy 
models. Indeed, dark energy models with a runaway potential of the Ratra-Peebles type 
for a real scalar field x suffer from the extreme smallness of the mass of x : today m x = 
O(Hq) ~ 10~ 43 GeV [6]. Such a small mass would leads to O(l) deviations from Newton's 
law over scales smaller than hc/m x ~ O(Gpc), if the scalar field x couples to matter with a 
strength similar to that of gravity. Local tests of gravity strongly rule out such deviations 
over scales larger than 0.1mm [1]. Generally these constraints are taken to imply that if x 
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does couple to baryonic matter it does so with a strength much less than gravity. This result 
is valid when the model does not have a chameleon property whereby the mass of the scalar 
field becomes environment dependent [7-10]. In particular, chameleon models evade solar 
system tests thanks to a thin shell effect with a large suppression of the scalar field force 
for large enough bodies such as the Sun. Of course, it could also be that baryonic matter 
does not couple to dark energy at all, which would still allow for non-trivial interactions 
between cold dark matter and dark energy [11—17]. In this case, gravitational experiments 
in the solar system would not be influenced by the scalar field at all. In fact, the thin shell 
effect of chameleonic theories implies that this is effectively what happens for structures up 
to galaxy clusters. In a sense, the coupling of the scalar field has been effaced by the thin 
shell property and an effective model emerges where large scale fluctuations in dark matter 
are the only 'species' which couple to dark energy. Gravity could also be modified on very 
large scales as in the DGP model on the accelerating branch (see e.g. [18-20] and [21] for 
a review). In this case, acceleration is entirely due to the modified structure of gravity. 

At the background cosmological level, models of modified gravity such as f(R) theories 
and its chameleonic siblings are extremely close to a ACDM model. This is also true of 
large scale modified models of gravity such as DGP. Testing the validity of these different 
approaches and distinguishing them can be envisaged at the perturbative level (for an 
investigation about the growth of perturbations in f(R) models, see e.g. [22,23]). Indeed, 
structure formation is sensitive to fine details in the modification of gravity. In particular, 
the Newton potential in these models does not remain unique but is embodied in two 
realizations ip and 4> with very different roles. First of all, the propagation of light signals 
is affected by ip + (f>. This is important when considering weak lensing or the ISW effect. 
On the contrary, peculiar velocities are sensitive to ip while the baryonic density contrast 
is related to 4> via the Poisson equation. All in all, depending on which correlation between 
lensing experiments and either peculiar velocity data or galaxy counts has been analysed, 
one may find different guises of the modification of gravity. Due to the coupling of the scalar 
field to matter, the growth of the CDM and baryonic density contrasts is also modified. 

Constraints on the coupling of baryons and dark matter to dark energy have been 
imposed on various scales. The tightest results have been obtained using CMB data for 
models with a coupling to dark matter only (for a recent analysis, see [24] and also [25,26] 
for a discussion of coupling dark energy and neutrinos to dark energy). It has been shown 
that the deviation of the effective Newton constant to the baryonic one cannot exceed 5 
percent at scales of order 10 Mpc while the ratio G e s/G could be as large as 2.7 for scales 
of 1 Mpc. A tighter bound coming from tidal disruption of dwarf galaxies amounting to 
a deviation of less than 4 percent at scales of order 100 kpc has been quoted although 
a comparison with data needs to be performed [27, 28]. The modification of Newton's 
constant for models with a coupling (3dm to dark matter only reads [25,30] 

where oldm = ^Pdm-i ^ ^ s ^ ne co- moving wave number and A c the co-moving range of the 
scalar field force which we will call the Compton length in the following. For scales larger 
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than the Compton length, we have G e g ~ G while well inside the Compton scale we find 
G e s « (1 + «bm)G. CMB data are therefore compatible with a Compton length in between 
1 Mpc and 10 Mpc. 

Deviations from General Relativity can be parametrized in a model independent fash- 
ion in terms of the slip function E and the growth rate, 7, (see e.g. [29]). These are defined 
by: 

V 2 (0 + V) = 3Zn m H 2 5, (1.2) 
In/ din 5 

t = r~Fr> f = Ti — L3 ) 

In i l m a in a 

where S is the density contrast. The definition of 7 is equivalent to / = ClZi- Alternatively 
one could work in terms of r/ rather than E where 77 = (f)/ip. These two definitions are 
equivalent when E and 77 are constant and 2E = 1 + rf x . In General Relativity, E = 77 = 1 
and 7 f» 0.55. 

In realistic models the Compton length varies in time, as the (effective) mass of the 
scalar field evolves in time. In some models such as the chameleonic ones, the Compton 
length decreases with time. It may happen that galactic scales enter the Compton length 
at a redshift z* and therefore the growth of structures for such scales is subsequently 
affected. In these models E and 7 would deviate from their GR values only at late-times 
(z < z*). Moreover, these parameters become time-dependent. Hence our analysis of 
such late-time modifications of gravity in scalar-tensor theories provide a template for the 
redshift dependence of E and 7. This could be useful for future observational studies of 
modified gravity. 

In the following, we will analyse the growth of structure in linearly coupled models when 
scales enter the Compton length and gravity is henceforth modified. The paper is organized 
as follows: In Section 2 we discuss the perturbation equations and the evolution of density 
perturbations. In Section 3 we discuss what happens if a cosmological perturbation crosses 
the Compton length. We study the growth of perturbations across the jump and find 
analytical expressions for the perturbation growth. The phenomenological consequences 
are explained in Section 4. Our conclusions can be found in Section 5. 

2. Cosmological Perturbations in Linearly Coupled Scalar Field Models 
2.1 Cosmological perturbations 

We are interested in the growth of cosmological perturbations in models where matter 
interacts with a scalar field. For convenience we work in the Einstein frame where the 
gravity equations take the usual form. This scalar field can mediate a new force between 
matter species which is non-universal. We are working in the Newtonian conformal gauge 
in which the metric reads: 

ds 2 = -a 2 (r)(l + 2^(x))dr 2 + a 2 (r)(l - 20(x))dx 2 . (2.1) 

The Einstein equation 

R»u ~ = SirGT^ (2.2) 



- 3 - 



has the usual form and the total energy-momentum tensor is: 

T,u = ^ 4 " AX P(aA A H A) (2-3) 
A 

+v M xV,x - \g^u [(v x ) 2 + m 2 X 2 ] . 

where A labels the different matter species and u^u^g^ = —1. We have also introduced 
K4 = v87rG. Notice the explicit coupling of the scalar field to matter. The couplings /3a 
can be all different from each other. We assume that the different matter species interact 

only gravitationally. Matter is conserved implying that V ' ^{p(A) u (A)p) = 0- Taking the 

{A) 

divergence of T^ v and requiring it to vanish in order to satisfy the Bianchi identities gives 
the x field equation: 

V 2 x = m 2 X + ^AP { A)^ AX . (2.4) 
A 

The acceleration of test particles ajf^ = u^ u Vuujf is influenced by the presence of the 
scalar field and reads 

a W = - K ,f3 A D^x, (2.5) 

where = V M + u^u" A ^ u V v . We also define for convenience piA) = eKi ^ AX P(A)- Since 
we are in the Einstein frame, the matter density p^) is the Einstein frame matter density 
which is not generally conserved. 

We define u^ 1 = v\ A Ja and then take Uq > and assume (v^l 2 -C 1 and tp, <ft <C 1 
and e <C f so that to leading order in each term: 

u^wa^l-^ + ^.Bj/ (2.6) 
Hence, including all potentially leading order terms, we have the Euler equation from Eq. 

cT l v\ A) + Hv\ A) + a-^djvl^ « -a' 1 [V + A x],i ■ (2-7) 

where H = a /a 2 . Notice that matter feels ip + /3aX an d not the Newton potential ip only. 

We also define pa = p A e SA /a 3 , and p A = const. Then keeping all potentially leading 
order terms, the conservation equation reads: 

oT 1 ^ + a- l v\ A) 5 A ,i » -a~ x div\ A) (2.8) 

We define T~L = aH and then consider perturbations to linear order. Using the divergence 
of the velocity field 6a = v l , A \ i we have: 

Sa ~ -0 A - (2.9) 
Oa + Wa ~ -bP + PAX]u = -i>AM. (2.10) 
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Now to leading order on sub-horizon scales we have, linearizing the matter perturbations: 

<f>,u ~ 47rGa 2 ^ P(A)$A, (2.11) 
A 

(2.12) 

X,u ~ a 2 m 2 x + a 2 ^ K4/3aP(A)^a (2.13) 

where the last one holds because we assume the energy density of the scalar field pertur- 
bation to be small. We have therefore: 

h + ns A ~ i)A,u - W - mj>. (2.14) 

On sub-horizon scales the last-terms can be ignored. Performing a Fourier transform we 
arrive at [25,30,31]: 

6 A + Wb A « -k 2 ^ A (2.15) 
= -U 2S ^^Ib{o)8b{1 + a A B{arn/k)), 

B 

where 

a AB (x) = TT ^. 

and x = am/k. The modification of gravity is sensitive to the Compton wave-length 
A = 1/m. Typically when the co-moving wave-number k becomes larger than the inverse 
Compton length A" 1 = ma, the growth of structures is affected [32]. 

Typical models with this behaviour have been constructed. Originally, the mass vary- 
ing dark matter model of [13] is such that only dark matter has a coupling to x- Moreover, 
the mass m increases with time (see also [33] for a model with growing neutrino mass). As 
a result for a given co-moving scale L = k^ 1 , gravity would have been modified in the past 
and come back to normal once L <C A c . This behaviour is also characteristic of models 
with a constant mass m. On the other hand, chameleon models [7,9,31] are such that all 
species couple to the scalar field at the linear level. At the non-linear level, a thin shell 
appears which suppresses the coupling of x to baryons. Phenomenologically, this can be 
implemented in this setting by imposing that only cold dark matter couples to x f° r scales 
below the size of clusters. The Compton length can increase or decrease in time. Models 
where it decreases can be obtained with an inverse power law potential of index n for which 
ma ~ a~( ra + 4 )/( ra + 2 ) decreases with time. In this case, gravity is modified at late time when 
scales become within the Compton length. 

2.2 Growth of structures: 7^ and jb 

In this section we focus on a two species system A and B, and £Ib 3> ^A so the matter 
density is then f2 m ~ $7b- In practise B would be Cold Dark Matter (CDM) while A would 
be the baryons. The evolution of structures is governed by 

Sa + U5 A w ^H 2 n B (a)5 B (l + a AB ), (2.16) 

5 B + Uh B ~ ^n 2 n B (a)5 B (l + a BB ). (2.17) 
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We define a(x) = (xbb(x) and 

1 + a AB (x) 



1 + Q-bb{x) ' 



When x > 1 or i < 1, ( is a constant. Hence whenever the scales are either well within, 
or well beyond the Compton length, £ ~ const and it is clear then that 

5a = -^-+ const, (2.18) 
1 + 4; 

satisfies the same evolution equation as 5 b- Provided the scales of interest do not cross the 
Compton length (so £ « const), at late-times we can assume 5a = 5b so 5a = (1+£)#b+Ao; 
Ao is a constant. At late times then 5a ~ (1 + £)#b- More generally we can define an 
effective value of £ thus: 

Sa = (1 + &ff)«*fl + A . (2.19) 
It follows that the bias, bAB = 5a/ 5b, between species A is given by: 

b A B = (l + £eff) + ^. (2.20) 



At late-times Aq/5b —> and so the dominant contribution comes from 1 + £ e ff • We note 
that if olab = olbb then £ = and £ e g = 0. In this case, 5a = 5b + Ao which is precisely 
the same equation as one encounters in GR. 

We suppose that the background cosmology is well approximated by the ACDM evolu- 
tion; in chameleon models, this is realized when m? H 2 . This implies that /5b, pa oc a -3 
and furthermore since £l m Ob S> &a'- 

n B ^n m = — 3 1 , (2.2i) 

where r2 m o is the value of Q m o at a = 1. 
We define: 

f B (\na,k) = ^^-, (2.22) 
dlna 



and find that this satisfies 



- ^f^l fB + fl + f B ~ |0B(a)(l + qbb). (2-23) 
where ' = d/dhxa. Henceforth Oa does not enter the field equation and we replace Ob by 

r2 m . 

We now assume that we know the solution, /gr, of Eq. ( |2.23| ) when «bb = i.e. in 
General Relativity. We parametrize the solution in the modified gravity scenarios (ctBB 
0) thus: /(lna, k) = (1 + <7B(lna))/GR,(lna). This gives: 

/ ~i I \ . 3f2 m , x 3aBB^m ,„ Q/A 

(1 + 9B)9B + TTP^dB + g B f G R = -^72 ' ( 2 " 24 ) 

z /gr z ^gr 
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In the ACDM model it was found that / GR « n° m 55 [35]. Since U m « fi B , ^b// G r ~ ^m 0,1 , 
which varies only very slowly up to the present epoch (a different approximation was be 
found in [36], but for the purpose of this paper, we use the one found in [35]). We find the 
analytical solution to Eq. ( 2.24 ) under the approximation ^m//cR ~ 1 i- e - /gr ~ ^m 5 - 



For the time being, we also approximate obb by a constant. We then have: 



5 25 3 , 

9B = S^Bb) ~ -J + y — + -«BB- (2.25) 

Thus with /gr = n2? R ~ ' we have / = Q]^ where: 

( i,\ ~ ln(l + g B (a,k)) 

7b(o, fc) - 7GR ~ r-p: • (2.26 

in il m 

We could similarly define gA and ^a to describe the evolution of the ^4-type matter species. 
Under the assumption that couplings are constant, we have 5a = (1 + £)<5b + Ao and hence: 



(1+5A) = (1 + 9b) 



bAB 



(2.27) 



At late-times the term in square brackets tends to 1 so gA — > 5b- 

However, as we shall see, if the Compton wavelength is crossed and £ ^ (i.e /3a ^ Pb) 
then the simple correspondence of Eq. ( 2.27| ) is broken and it is entirely feasible that 75 



could deviate greatly from 7gr whilst 7^4 hardly changes at all. 

In general though, we would have 7a, 7_B 7^ 7gr if «bb 7^ 0. We study the effect of 
the jump in obb an d £ when scales get inside the Compton length in the following section. 

2.3 The Slip Functions 

Being interested in the modification of gravity, we now focus on the two gravitational 
potentials, </> and ip and their ratio described by the r\ = <f>/ip parameter or 2S = 1 + g^ 1 . 

The quantity ip + <f) is invariant under a conformal rescaling of the metric, however, 
individually, <f> and tp are not. A choice of conformal frame is essentially a choice of standard 
ruler and clock. In General Relativity there is a preferred and obvious choice of conformal 
frame where the Newtonian constant, G, is fixed and the energy momentum tensor is 
conserved ensuring that the masses of particles are constant. In modified gravity such a 
frame choice is not generally possible. In some theories (i.e. those with a universal coupling 
to a scalar field), there exists what is commonly known as the Jordan frame. In this frame 
the matter energy momentum tensor is conserved, particle masses are constant and non- 
gravitational physics is independent of space-time position. It is for this reason that the 
Jordan frame is often referred to as the "physical frame", however this nomenclature can 
be misleading. Additionally in the Jordan frame the effective Newton constant varies with 
space and time. In scalar-tensor theories, including those with multiple scalar fields and 
different couplings, one may always define an Einstein frame where G is constant, but 
T£;jj 7^ 0. In this frame the particle masses depend on space time positions and hence so 
does local non-gravitational physics. However, this does not mean that the Einstein frame 
is in some sense 'unphysicak. Both Jordan and Einstein frames, and indeed any other 
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choice of conformal frame, are physical in the sense that provided one does not assume 
or require quantities to be constant that are not and interprets all measured quantities 
correctly, they represent a perfectly accurate descriptions of nature. 

2.3.1 Weak Lensing Measurements 

We first consider weak gravitational lensing measurements, which are sensitive to (j) + ip 
(for a discussion of weak lensing in scalar tensor theories, see e.g. [34]). In the Einstein 
frame V 2 (j) = iirG5p, if one assumes that this is generally true and takes ip = r\~^(p where 
r] is assumed to be constant, we have: 

V 2 (<p + ip) = 8ttGZ5 Pi (2.28) 

where 2X = 1 + r] . In modified gravity models such as the ones we are considering, 
the parameter X defined by (2.28) is not constant. In general, we will take this form of 



the Poisson equation in the Einstein frame as the definition of X. In General Relativity, 
we have X = 1. In the type of scalar-tensor theory that we are considering, the absence 
of anisotropic stress in the Einstein frame implies that X = 1. Let us consider another 
frame obtained by a Weyl rescaling of the Einstein metric; g^ v — > g^ u = e~ 2w g^ u where 
w <C 1 is an arbitrary function. Under this change of metric, the Newton potentials are 
transformed as ip — > ip + w and <j> —> <j> — w, implying that <p + ip is frame invariant. This is 
not the case of the energy density of non relativistic matter p = —g^T^ which transforms 
as p — > p = e Aw p while G = e~ 2w G . Using V 2 = e~ 2w V 2 and the fact that for a given 
function w which is not a field itself bp = e~ 4w 5p, we find that X is invariant under a Weyl 
rescaling of the metric. This is a major advantage of defining X using the modified Poisson 
equation. Notice that when the Weyl transformation is designed to efface the coupling of 
a particular species A to x, on e must choose a field dependent 

w( X ) = -kaPaX, (2.29) 

as it is the case for the theories considered here. The resulting metric g^ v is the Jordan 
metric for species A. Unless the couplings are universal, this frame is not the Jordan frame 
for species B ^ A. In this frame we have 5p = e^ K ^ AX (bp + \k^a^XP)- Hence X is frame 
independent in this case if and only if K4/3a|<5x| ^ " ^ ,e - ^ * ne sca l ar field fluctuations 
can be neglected. When the coupling is universal and upon using ( 2,13j ), this is always 



true as long as /3 2 £1 C dmH 2 < m 2 ^i . This gives a bound on f5 depending on the ratio 
m/H S> 1. When this bound is satisfied, the slip parameter is frame independent. 

We denote X defined by Eq. (|3|) by X K . Provided one can measure 5p, X K can be 
extracted from weak-lensing measurements. Fourier-transforming Eq. ( [2.28 ) gives: 



-k 2 ^ + ijj) = 3^ 2 X K n c5c ~ m 2 n m Z K 5 B , (2.30) 



c 



3U 2 n m [E«6^] 5 A . 



where we have used the definition of the bias function in the last equality. Using weak 

S 2 
'b 



lensing measurements alone one can compare the ratio of X 2 <5o with its value at an earlier 
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time, which provides a measurement of the growth rate of T, K 8 B . We define 



S B (z,k)=D B {z,k)5i(k), (2.31) 

so din Db /din a = f B and <5, is the primordial value of the perturbation. The initial 8f(k) 
is proportional to the primordial power spectrum and is measured by CMB experiments 
such as WMAP. Therefore by combining weak-lensing measurements with an ansatz or 
measurement of the primordial power-spectrum one measures not T, K (z, k) but the combi- 
nation: 

X KK = X K (z,k)D B (z,k). (2.32) 

Another method of extracting S K would be to directly measure 8 a- For instance if A 
represents galaxies this could be done using galaxy surveys. From cross-correlation of 
weak-lensing and 8a one can then extract the quantity: 

H KA = Z K (z,k)b AB (z,k) (2.33) 

Combining measurements of 5a with a measurement of the primordial power spectrum 
provides: 

D A (z, k) = b AB (z, k)D B (z, k). (2.34) 

We note that D B and Da are determined by j B and ja respectively and by the cosmological 
model through U m (a). If one assumes a GR growth rate for species B then the measured 
value of S K is S Km : 

D B 

S Km = ^st; — • (2.35) 
^GR 

Alternatively if one assumes a particular ansatz for bA B = b(z, k) then one measures: 

S k6 = E K -^-. (2.36) 
bAB 

In a scalar-tensor theory, such as the class of theory considered here, £ K = 1. However 
since one does not measure S K on its own, the measured value of S K will depend on the 
ansatz one makes for either the growth rate 7b or bias bAB as well as the cosmological 
model, for instance: 

y .2*. y - 1 

Z-'Km — j ^nb 



£>GR &AB 

If one estimates b by fitting to the observables and assuming a GR growth rate Da = 
b AB D B = bD GR would give: b/b GR = D B /D GR and S Km = 

Given that weak-lensing measurements directly probe ip + (f>, one could measure r\ = 
<j)/ip directly if one could measure either (f> or ip directly. In this case, the frame in which one 
effectively measures r] depends on the method one uses to measure <fi or ip. The most direct 
method for measuring ip, which does not depend on the bias between different species, is to 
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use peculiar velocities. From Eqs. ( |2.10| ) and ( 2.1 3| ) it is clear that the peculiar velocities 
of species A depend on 

1p A ~ "(I + «AB) 3 ^ $B = ^(1 + Ol A b) [<t> + i>] ■ 

If /3a = @b so uab = olbb = ct, then ipA is precisely the value of ip one would calculate in 
the Jordan frame. Similarly if (3a = 0, then tpA is equal to the Einstein frame value of ip. 
Assuming that if) = ipA one would therefore estimate rj = r]g where: 

1 + V6 = —, — = t— • (2.37) 

ipA 1 + aAB 

We may also define 2T,g = l + r/g 1 = 1/(1 — «ab). In any scenario where species A does not 
couple to the scalar field, /3a = 0, T,q = t]q = 1. In theories such as f(R) theories, galaxies 
are expected to be effectively decoupled from the scalar force as the result of the chameleon 
mechanism which allows for compatibility of the theory with local tests of gravity. In such 
theories then measurements of the peculiar velocities of galaxies would find £ = r] = 1 and 
hence not reveal any modification of gravity even though large-scale CDM perturbations 
might well feel an non-negligible fifth force. 

2.3.2 ISW Measurements 

The ISW effect depends on the quantity <f)' + ip' . Given S K as defined Eq. ( |2.30|) is follows 
that: 

-k 2 [4/ + « 3n m n 2 ^! [/gr - 1] S B = 3n m ft 2 £ 7 &^ [f GR - l] 8 A , (2.38) 

where 



i-/b + K 



1 - /gr S k 



(2.39) 



We note that in the simple case where fs = /gr an d £ K = const, £/ = S K . As with 
weak-lensing measurements, in practice one would not measure £7- directly but either: 

£77 = ^iDb, 

or 

^IA = ^ib~AB- 

If one assumes a GR growth rate, then the measured of £/ would be: 

Db 

£/m = £/— -. (2.40) 
^GR 

In the class of scalar-tensor theories we have considered £ K = 1 and so: 

v (1 - /b)£>b (0 ... 

E/m = (1 - /gr)^gr (2 - 41) 
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2.3.3 Summary 

We have introduced the slip parameter E K using a modified version of the Poisson equation. 
Extracting the slip function from data requires an assumption on the growth factor. Using 
weak-lensing and assuming a GR growth rate one would measure: 

Db 



'-•Km — t-. ; 
^GR 

Using ISW measurements, and under the same assumptions one would extract: 

(1 - fn)D B 



J Im 



(1 - /gr)-C>gr 



Both expressions are determined by D-q which itself is given by /b = (1 + 5s)/gr and 
the growth rate 7b- The bias between species A and B depends on £ e g- and gs as does 
the growth rate of species A; ja- In these next section we will find very good analytical 
approximations to the evolution of gs and £ e g when the couplings /3a. and /3b jump at 



3. Effect of a Change in the Coupling 

We now consider what happens when the values of cxab and qbb change at some redshift 
z = z* for instance between a perturbation crosses the Compton wavelength of the scalar 
field x- We first consider the growth function for the dominant matter species B before 
considering the sub-dominant species A. 



3.1 Dominant Species Growth Function, fs 

To simplify the problem, we assume that the coupling function obb changes abruptly 
(jumps) at a redshift z*, i.e. we assume that for z > z*, obb = and for z < z*, 
otBB = oi for some a and «o- We define go = g{ao) and g a = g(a). For z > z*, o-bb is 
taken to be constant back to the far past, and so gs = go = const. For z < z*, gs must 
satisfy Eq. ( gjg ): 



, \ ^^^772 9 Ft 

k 1 + gB)gB + ttp)-9b + " 



V? 



GK 



/gr 



J 



GR 



(3.1) 



We define gs = g a + Ag and using Eqs. and ( gjg ) arrive at: 



/gr 



3Q r 



2/« 



+ 1 + 2g a + Ag 



GR 



Ag. 



(3.2) 



f^ 55 . As in the derivation of g(a), we take Q//q R 



Now / GR = 

(and constant). This is equivalent to taking /gr 



£2 m to be unity 



/o = yHm for the purposes of finding 



gs, £ e ff and gA- Clearly this is exact when VL m = 1 and is a good approximation generally 
when O~ 01 — 1 <C 1. Notice that /gr ~ /o = &m is the only approximation that we make 
in what follows. 



- 11 - 



Growth rate, f , for (z t , a, aj = [1 , 1 , 0] Growth rate, f B , for (z + , a, aj = [2, 1, 0] 




Redshift: z Redshift: z 



Growth rate, f B , for (z + , a, aj = [4, 1 , 0] Growth rate, f B , for (z + , a, a Q ) = [1 , 6, 0] 




Redshift: z Redshift: z 



Figure 1: Sample behaviours for the linear growth rate, /s, of the dominant matter species i.e. 
dark matter. In all cases shown above the ao, fifth force coupling for z > z*, is taken to vanish 
and the late-time coupling a, for z < z», is a constant. The plots above show the behaviour for 
(z*,a) = (1,1) (top-left), (2,1) (top-right), (4,1) (bottom-left) and (1,6) (bottom- right). We see 
that in all cases there is a pronounced deviation from General Relativity where Jb = /gr ~ ^m 545 
(shown as the dot-dashed blue line). The solid red line is the exact numerical solution, and the 
dotted black line is our analytical approximation: f B = (1 + i g B )f2^ 545 . We note that the analytical 
approximation derived in this work is almost exact at all times. 



With this approximation we find that we can further set 3Q m /2fo+l + 2g a 5/2+2# Q 
is a constant and we define it to be hq. By integration we find: 



/x + Ag 

where Aq is a constant of integration and Dq is defined by 

D' = f D 

and at Dq(z^) = D*. We have 



D = D, 



K(n r , 



i 



1/3 



(3-3) 



In this equation we have defined = Q. m (z = z*). As f2 m ,f]^ — > 1, Dq/D* 
(l + z*)/(l + z). 



(3.4) 
a/a* = 
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Now at z = z*, gs = go and so Ag = (go — g a ) and, hence 

9a - go 



A 



vo - {g a - go) 



Thus we find 



9B ~ g a 



/j. A D 



Dn. + A D%° ' 



(3.5) 



(3.6) 



9b ~ 5o + (5« ~ go)G B {Do/D*;Vo,A ), 
where for X < 1, Gb(X; /x , Aq) = and for X > 1: 

" X' 10 - 1 " 



G B (X > 1;mo,A)) 



IN + At 



Now: 



[lnD B ]' = (l+5B)[ln£> GR ]'. 
Thus it follows that using Z?gr ~ -Do that: 

« D 9 °F B (D /D,;fi ,A ) 



where we assume that at some initial time, z 
X < 1, F B (X; no, Aq) = 1. For X > 1. 

F B (X>l;/xo,A )=X^ 



(5,- and D ( 



o 



£>b 



-do 



1 + iol-w 
1 + 4) 

Excellent approximations to both <5b and /b are then given by: 

/b = (1 + 5b)/gr = (1 + 5 b)^ gr " - 55 . 



(3.7) 
(3.8) 

1 while for 
(3.9) 



(3.10) 
(3.11) 



These analytical approximations of both /b and <5b become exact in the limit f2 m — >• 1. 

Using this analytical approximation we can straightforwardly calculate the growth 
index: 

ln(l + <? B ) 



IB ~ 7GR + 



ln$7n 



(3.12) 



however since /b does not tend to 1 as Q m — > 1, we follow Ref. [37] in noting that it is 
arguably better to parametrize /b as (1 + <7B)^n? than as fim . In the latter case 7b often 
diverges as fi m — > 1 but /b remains finite. In the former case 7b = 7gr ~ 0.55 and gB is 
as we have calculated. 

Figure (|l|) compares of the exact value (solid red line) of fg, calculated by numerically 
integrating the perturbation equations, and the analytical approximation (dotted black 
line) presented above for given values of a, Qo and z*. We also show the GR growth rate 
as a dot-dashed blue line. We see that in all cases our analytical approximation provides 
an excellent fit to simulations. We note that because a > 0, fs > /gr- 
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3.2 The Slip Functions: T, Km and T,j m 

We noted that if one measures the slip function, S, from either weak-lensing or ISW mea- 
surements assuming that the growth rate is unaltered then what is measured is respectively: 
S rem and E/ m . In these theories S Km = Db/Dqr, and the approximate analytic form of 
this for a coupling that changes at z = z* was given above in Eq. (pT 



Weak Lensing Slip Parameter, £ , for (z, a, a ) = [1 , 1 , 0] 



Weak Lensing Slip Parameter, £ , for (z t , a, a ) = [2, 1 , 0] 





Numerical solution 

Analytical Solution: Z K = D B /D GR " 
GR slip parameter (=1) 









Redshift: z 

Weak Lensing Slip Parameter, £ , for (z , a, a ) = [4, 1 , 0] 



Numerical solution 
Analytical Solution: I 
GR slip parameter (=1) 
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Numerical solution 

Analytical Solution: E K = Dg/D GR 

GR slip parameter (=1 ) 









Redshift: z 

Weak Lensing Slip Parameter, £ , for (z , a, a ) = [1,6, 0] 




- Numerical solution 
Analytical Solution: I 
GR slip parameter (=1 ) 



Figure 2: Sample behaviours for the weak-lensing slip parameter, E Km = Db/Dgr, derived under 
the assumption of a GR growth rate for perturbations. In all cases shown above the ao, fifth force 
coupling for z > z*, is taken to vanish and the late-time coupling a, for z < z*, is a constant. The 
plots above show the behaviour for (z*, a) — (1, 1) (top-left), (2, 1) (top-right), (4, 1) (bottom-left) 
and (1,6) (bottom-right). In GR, S Km = 1 at all times (shown as the dot-dashed blue line). We 
see that in all cases, S Km grows monotonically from an initial value of 1 for z < z*. The solid red 
line is the exact numerical solution, and the dotted black line is our analytical approximation. We 
note that the analytical approximation derived in this work is, as with the approximation for the 
growth-rate, almost exact at all times. 

We have also calculated all that is required to have an analytical approximation of the 
measured ISW slip function S/ m since S/ m /S Km = (1-/ b )/(1-/gr) and / B = (1+5b)/gr 
so: 

S/ m /S Km = 1 - f B _ ~ 1 - O _ 5 55 _ -i • ( 3 - 13 ) 
J GR L 12m L 

Now for qbb > 0, <7b > 0, and since /gr < 1 for f) m < 1 it follows that E/ m < S Km and it 
is possible to have S Km > 1 and S/ m < 1 at the same time. Additionally, as £l m — > 1, S/ m 
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can grow large and negative unless g-g decreases at a sufficient rate. It is feasible that X 
would show little deviation from its GR value whilst |S/ m — 1| ~ 0(1)- 




012345 012345 
Redshift: z Redshift: z 



Figure 3: Sample behaviours for the ISW slip parameter, £/ m , derived under the assumption 
of a GR growth rate for perturbations. In all cases shown above the ao, fifth force coupling for 
z > z*, is taken to vanish and the late-time coupling a, for z < z*, is a constant. The plots above 
show the behaviour for = (1,1) (top-left), (2,1) (top-right), (4,1) (bottom-left) and (1,6) 

(bottom- right). In GR, E Km = 1 at all times (shown as the dot-dashed blue line). The solid red 
line is the exact numerical solution, and the dotted black line is our analytical approximation. We 
note that the analytical approximation derived in this work is, as with the approximation for the 
growth-rate, almost exact at all times. 

Figures (|2|) and @ respectively compare of the exact value (solid red line) of S Km and 
S/ m , calculated by numerically integrating the perturbation equations, with the analytical 
approximation (dotted black line) presented above for given values of a, ao and z*. We 
also show the GR value of S Km = S/ m = 1. In all cases it is clear that our analytical 
approximation represents an excellent fit to simulations. We also note that for a > 0, 
S Km > 1 and grows monotonically for z < z*. S/ m , however, can be both > 1 and 
< 1, although typically > 1 values are only found when is sufficiently in the past. 
Additionally E/ m can become negative. Typically the deviation of E/ m from 1 is noticeably 
more pronounced than that of 
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3.3 Fifth Force Linear Bias: b x 

We define the quantity b x so that 

5 a = b x 5s + A = bAB^B, 
where Ao is a constant and bAB is the linear bias. We define the usual linear bias: 

A " 



hm(SA) 



1 



In the absence of a fifth force (i.e. in GR) we then have simply bAB = b\m- in the presence 
of a fifth force, we have instead 

b A B(k, z) = b x (k, z)bi in (5 A (k, z)), 

and so we may see b x as being the additional contribution to the linear bias due to the 
fifth force. 

If one measures 5a then using the form of bu n and the fact that b^(5A)8A is pro- 
portional to the initial Gaussian perturbation 5i, both in GR and this modified gravity 
scenario, one can measure Ao and hence by m directly using the bispectrum of 5a per- 
turbations. It appears as a shift in the bispectrum of 5a compared to the primordial 
spectrum. It is therefore possible to measure the 'bias-corrected' value of 5a which is given 
by 5a^ = ^lTn^ = bx^B- Thus bias due to b x remains even when one corrects for the usual 
linear bias. 

If species A represents galaxies there will generally be additional multiplicative con- 
tributions to the bias relative to the linear CDM perturbation 5b due to, for instance, the 
rate of galaxy formation and non-linear effects. In many viable modified gravity theories, 
modifications to gravity are suppressed on short-scales and when the ambient density is 
much larger than the cosmological density. This is the case for chameleon models. So 
we estimate that the formation and non-linear bias would not be greatly affected by the 
presence of a cluster scale modification of gravity. Hence b x is therefore expected, in many 
theories, to represent the dominant additional contribution to the bias due to fifth force 
effects. 

If «ab = «bb at all times (even if «bb / 0) then b x = 1. A non- unity value of b x only 
results when species A couples to the fifth force with a different strength than does species 
B. 

We now calculate b x across a jump at z = z*. For z > z* we define: 

b x o = (1 + cxab(z > z*))/(l + &bb{z > z*)), asB = ao, 

and for z < z*, 

b xf = I 1 + ®ab{z < 2*))/(l + oibb{z < 2*)), a B B = a. 

For z > z*, we have simply b x = 6 x o- I n the far future of the transition at z = z*, b x — > b x f 
but generally b x (z < z*) / b x f. 
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Fifth Force Linear Bias, b , for 6 =0 and (z , cc, a ) = [1 , 1 , 0] 



Fifth Force Linear Bias, b , for (5 = and (z , a, a ) = [2, 1 , 0] 
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Analytical solution 
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Figure 4: Sample behaviours for the fifth force linear bias parameter, b x — S A /8b- The usual 
linear bias bn n depends only on 5a and can be measured directly using higher statistics of 6 a- In 
all cases shown above an, the fifth force coupling for z > z*, is taken to vanish and the late-time 
coupling a, for z < z„, is a constant. The plots above show the behaviour for (z», a) = (1, 1) (top- 
left), (2,1) (top-right), (4,1) (bottom-left) and (1,6) (bottom-right). In GR, b x = 1 at all times 
(shown as the dot-dashed blue line). We see that in all cases, b x decreases monotonically from an 
initial value of 1 for z < z* . The solid red line is the exact numerical solution, and the dotted black 
line is our analytical approximation. We note that the analytical approximation derived in this 
work is almost exact at all times. 



For z < z», we define Aft = b x — b x f and then inserting 8a 
evolution equation for 5a we arrive at 



Y 



3fl r 



Y' 



+ An in to the 



(3.14) 



where Y = {M))8b/8* b . We define S = Y' / JgrY and using the above identity we have 



/gr 



3Q r 



i + s 



GR 



s. 



(3.15) 



As with the deviation of the g-Q expression, the only approximation we make is to assume 
/gr ~ fo = ^m 1/2 - Then we 3Q, m /2f$ -1 = 1/2 and have: 

S' _ D'o 
a + S)S Do' 
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and so: 



S 



Y' 



A, 



f Y 2(( J D / J D*) 1 / 2 -^i)' 
Continuity at z = implies = 1 + go and so 

1 + 50 



(3.16) 



A x 



Thus 



y « (6 x0 - 6 X /) 



where we have used y* = (6 x o — b xf)- Using the expression for 
fiQ = 5/2 + 2^0, then we obtain 



b x = ho + ( b x.f - ho) [1 - B ( D o/ D *;9o,9a)} , 



(3.17) 

derived above and 
(3.18) 



where 



BpC < l;go,5a) = 1, 

B(X > l;go,g Q ) = X^ 1+ ^ 



1 + A 



o 



1 + ,4o^ _/to 



1 _ AxX- 1 ' 2 
I -Ax 



(3.19) 
(3.20) 

This analytical approximation to b x is, as with the 5b approximation, exact when Vt m = 1 
and remains an excellent approximation up to the present era. Figure (Q) compares of the 
exact value (solid red line) of b x , found via numerically integration of the perturbation 
equations, with the analytical approximation (dotted black line) presented above for given 
values of a, ao and z*. We also show the GR value of b x = 1 as a dot-dashed blue line. In all 
cases it is clear that our analytical approximation represents an excellent fit to simulations 
at all times. For a > 0, b x < 1 and decreases monotonically for z < z*. 

3.4 Growth Function, /a, for Sub-Dominant Matter Species 

It may be the case that one tracks the growth of a perturbation in the dominant matter 
species, B, using a a sub-dominant matter species, A. It would therefore be important to 
know how the growth function of 5a relates to that of 5b- 

Now we have define the 'bias-corrected' value of 5a to be 5^ = ^^(^a)^; b im(^A) 
can be estimated directly from 5a measurements using higher order statistics. Thus whilst 
estimating bias directly from the data allows one to correct for bn n , the b x contribution to 
the bias remains. We therefore consider 



Abe) 



bJk,z)5 



B, 



and define: 



din ,5 



(be) 



din a 



(3.21) 



(3.22) 
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Growth rate, for (z , a, a ) = [1 , 1 , 0] Growth rate, f^', for (z #1 a, a ) = [2, 1 , 0] 




Redshift: z Redshift: z 



Growth rate, f, for (z t , a, og = [4, 1 , 0] Growth rate, f, for (z t , a, aj = [1, 6, 0] 




Redshift: z Redshift: z 



Figure 5: Sample behaviours for the linear, bias-corrected growth rate, f A , of a sub-dominant 
matter species which is uncoupled to the fifth force e.g. galaxies in chameleon / f(R) models. In all 
cases shown above, ao, the CDM fifth force coupling for z > z*, is taken to vanish and the late-time 
CDM coupling a, for z < z*, is a constant. The plots above show the behaviour for (z*, a) = (1, 1) 
(top-left), (2,1) (top-right), (4,1) (bottom-left) and (1,6) (bottom-right). We see that in all cases 
there is a pronounced deviation from General Relativity where f A bc ^ = /qr ~ ^m 545 (shown as the 
dot-dashed blue line). The solid red line is the exact numerical solution, and the dotted black line 
is our analytical approximation: f A = (1 + .g^f^ 545 . We note that the analytical approximation 
derived in this work is almost exact at all times. 



Thus it follows directly from Eq. ( |3.21 ) that: 



(be) 

9 A 



9B + 



b' 



(3.23) 



/grV 

For z > z*, gB = go, b x = b x o = const and so c '(z > z*) = go- For z < z* we note that 

(3.24) 



<J n 



&v = b xf + /y, 



where Y obeys Eq. ( |3.14| ) and so 



K ~ b Xf 



;i+5b) + 



Y' ' 

IgrY 



(3.25) 
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Making the single approximation /qr ~ f 
by Eq. ( [3.16 ) and so: 



/ f2 m we have Y' / fyY = S where S is given 



9a ( z > z *) =9o + {9a - go)G A (D /D*), 



where 



G A (X < 1) 
G A (X > 1) 



5 x/ 



I» - 1 



X»° + An 



+ 



5 X/ 



and /i 6c) 



1+ffo 

9a~ 90 



(1 + <m)/gr- We can also have = 6^ /Si in terms of Dq\ 



X 1 ' 2 - 1 
X 1 / 2 _ A x 



(3.26) 

(3.27) 
(3.28) 



(be) 



(3.29) 



Now when /3a = i.e. the A-type matter species does not feel the fifth force and we have 
l + £ = 1/(1 + q) and l + £ = 1/(1 + a ) and (£-£o)/(l + £) = -(a-a )/(l + a ). Thus, 
it can be checked that, approximately 



9 A 



9o 



9B - 9o 



l + £ \ x ? (l + 2x 5 ) 



l+£eff 



(3.30) 



where = (£ e ff — £o)/(£ — £o); this approximation is particularly good for small i.e. close 
to z = z* and at late times when xe ~ 1. We also find that (without requiring (3 A = 0), 
there is a similar approximation for gg — go in terms of xg. Specifically: 



gB- go - (g a -#o)tanh 



c x ? 



where 



3(a — ao) 



^2 ' 



(l + a)(g a -g ) 

Hence we can roughly relate g^' and hence ja to and hence b x : 

x$ (1 + 2x ? ) 



ffA ~ go 

9« - go 



(be) 

1+j 

1+Ceff 



tanh 



L — 



The analytical approximation of given by Eq. (|3.26j ) is exact at early times and remains 
very accurate to late times. Whichever approximation one uses, the growth function, 7^4, 
is then given by: 



1A ~ 7GR + 



ln(l + g A ) 
lnf2 m 



(3.31) 



All of the expressions derived have been found using only the approximation Q. m / ffi ~ 
f2~ 0,1 ~ 1 and hence share the property that when Q m = 1 they are exact. Fortunately 
they also remain excellent approximations up to the present day when O m ~ 0.26. Figure 
(P) shows for = at all times for the same values of a, qq and as were plotted 
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for fs in figure (|T|). In each plot the solid red line is the exact numerical solution, the 
dotted black line is the analytic approximation f A bc ^ = (1 + <7^)f2^ 545 with g A from 
Eq. ( 3.26| ) and the dash-dotted blue line is the GR growth rate, /gr- It is clear that our 
analytic solution represents an excellent approximation to the exact behaviour of f A c \ 

In all cases, we have considered ao = and the coupling to the dominant species (i.e. 
CDM) turns on at z = z*. If there were a universal coupling (Pa = Pb) to the fifth force 
/a°^ = /b as shown in Figure (|l]). By comparing figures (||) and (|5|), we note that when 
the coupling to the dominant B matter species turns on at late times, the deviation of f A c ^ 
from its General Relativity value is much less when species A is uncoupled /3a = than 
when Pa = Pb- We note that, whatever value Pa takes, a > 0, implies f A > /gr- 

We have considered a scenario where the dominant form of matter involved in large 
scale structure formation i.e. cold dark matter (CDM) on cluster, and larger, scales, is 
coupled to an additional fifth force mediated by a scalar field, x- We have denoted such 
matter to be species B and its coupling strength Pb to x- The fifth force between species B 
matter particles is then cxbb = ^Pb times the strength of gravity. One does not, however, 
observe CDM on cluster scales directly, only its effects on observable, non-dark, forms 
of matter. We have therefore also allowed for a second, subdominant, type of matter 
dubbed species A, whose density contrast 5a, and peculiar velocity, 9a, perturbations can 
be directly observed and are used to extrapolate information about the density perturbation 
of large scale distributions of CDM (i.e. species B). Typically, one uses the distribution 
and velocities of galaxies to measure the distribution of CDM, and so our species A should 
be taken to be the baryonic matter in the universe (i.e. galaxies). 

We have assumed that galaxies, as a whole, have a different (effective) coupling to % 
than CDM. We denote this coupling by Pa and defined olab = %PaPb- Since galaxies are 
also predominantly constituted of cold dark matter, in the simplest scalar field models (i.e. 
those with approximately linear field equations) Pa = Pb, and olab = olbb- However, 
there is a key difference between the CDM confined in galaxies and that distributed on 
large scales, namely the former is much denser than the latter and inside the galaxies 
the density perturbation compared with the background is highly non-linear. In scalar 
field theories with non- linear field equations, a chameleon mechanism might develop which 
causes the mass of the scalar field to depend on the environment. Typically x would then 
be much heavier in denser regions than it is in sparse regions. If the mass of the field, m x , 
inside galaxies is sufficiently large (i.e. rn^ 1 <C fewkpc), then the galaxies would effectively 
decouple from external perturbations in x an d so there would be almost no fifth force i.e. 
Pa <S Pb- This is realized in both chameleon / chameleonic f(R) theories (in which there 
is also a coupling to baryonic matter) and the related varying mass dark matter models. It 
is perfectly feasible therefore that one might find gab *C «bb at least when one considers 
the evolution of perturbations on scales where they are linear i.e. 5b <C 1. When the 
CDM perturbations go non-linear, it is feasible that they might also decouple from the 
fifth force in the same way as the galaxies have. This final possibility is beyond the scope 
of this work and, since it involves complicated non-linear behaviour, mostly likely requires 
TV-body simulations to address fully. 
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4. Constraints 

In this section we briefly describe the range of possible behaviours for the observables 7a 
and £ for different values of the couplings /3a and (3 b • We use recent constraints on the 
growth parameter and £ to give limits on the different couplings. Good limits come from 
measurements of the galaxy and Lyman-a (Lya) power spectra. These measure the ratio of 
the power spectra at one value of z to that at another. Neither of these tools directly probe 
the growth rate of the large scale cold dark matter perturbations. Lya absorption systems 
probe the power spectra of the density perturbations in baryons, and galaxy surveys probe 
the galaxy power spectrum. 

The growth rate determined from both galaxy and La power spectra are usually sub- 
ject to some bias correction. For the Lya systems this is generally done by comparing 
observations with iV-body simulations which assume GR. In galaxy surveys such as SDSS 
the bias is assumed to have a luminosity dependence which is specified a priori up to an 
overall normalization constant which is fitted for by comparing observations with the CMB 
linear power spectrum extrapolated to the current epoch (assuming GR growth). Neither 
method accounts for the additional red-shift dependent bias, b x , resulting from a non- 
universal coupling to the fifth force. We do, however, assume that any bias present in GR 
has been removed. Thus we take the measured density perturbation to be 5^ (up to an 
overall normalization constant which is degenerate with the average bias). The measured 
growth rate is therefore /j^ = (1 + <7^)/gr- 

Ref. [37] catalogues recent limits on from Lya systems and from the 2dFGRS 
galaxy survey. The 2dFGRS limit is f^\z = 0.15) = 0.49 ± 0.14 [38,39]. For this 
measurement the linear bias was estimated directly from the data essentially by assuming 
b^ 1 (SA)^A = D A bc \z)5i where Si is the initial Gaussian fluctuation. Thus the linear bias 
estimated in this case will be bo(8A) in both GR and the class of modified gravity model 
consider here, and the quoted value of / is truly f A b °\ 

Recently Ref. [38] reports another limit of this time at z = 0.77 of /bo = 
0.70 ± 0.26 (VVDS) from the VIMOS-VLT Deep Survey. This time, the linear bias, b 
was estimated using a different method. It was assumed that 6 cst = b^^ where 6 C st = 
a^ gal) {z = 0.77)/4 cmb) (z = 0.77) = 1.3 ± 0.1; a ( ^ mh)2 is the normalization of the power 
spectrum amplitude extrapolated from WMAP assuming a GR growth rate. Thus in this 
model b cst = b sa ^ / Dqr. The constraint on / quoted in Ref. [38] therefore corresponds 
to f^ ] D^ c) /D GR = 0.91 ± 0.36 at z = 0.77. 

The limits from Lya systems are all for z > 2, and the best constraints are S^ a (z = 
2.72)/4yi(z = 2.125) = 0.83 ± 0.11 from Ref. [41] and f Lya {z = 3) = 1.46 ± 0.29 [40]. 

In what follows we consider the constraints on /3a, 13b and .z* that follow from the 
limits on the galaxy / Lya perturbation growth rate collated. 

We consider the la constraints on Pa, /3b and z* that arise from limits on the deviation 
of and from their GR values. Here species B is the linear (large scale) cold dark 
matter perturbation, and A is either galaxies or Lya systems corrected for GR bias. 
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2dFGRS, Lya and VVDS Data 



Universal Coupling 


/3gal = 




a 




a 


0.5 


< 0.52 


0.5 


<57.4 


1 


<0.31 


1 


<6.7 


2 


<0.24 


2 


<1.6 


3 


<0.22 


3 


<0.88 


4 


<0.35 


4 


<1.3 



We begin by considering the case of a uni- 
versal coupling when Pa = Pb = P, that turns 
on at some redshift z = z*. So 2/3 2 (z > z*) = 
and 2/3 2 (z < z*) = a. We assume a flat ACDM 
background with f2 m = 0.27±0.1. The 1-u limits 
on a for different z* from the 2dFGRS, VVDS 
and Lya data are shown in left half of table [I]. 
These are derived by finding the a that mini- 
mizes the x 2 ( a ) f° r the data [37]; a say. The 
1-cj confidence limit corresponds to those value 
of a for which x 2 ( a ) ~ X 2 (&) < 1- For z* > 3, 
the constraints on a weaken slightly because of 
large value of /L ya at z = 3 extrapolated from 
the Lya data. 

Next we consider the limits on a = 2j3^{z < 
z*) under the assumption that galaxies effec- 
tively do not couple to the fifth force /3 ga i = 0. 
We assume that P^ya = Pb- Again for z > z#, 
Pb = and the coupling turns on at z = z*. 
Such a scenario would be realized, for instance, 

in f(R) and chameleon models where the effective coupling is density dependent and as 
such virialized objects such a galaxies would be uncoupled. These limits are given in table 
[l[ Table |2| shows limits on a for scenarios where the fifth forces couples only to dark 
matter and baryons are uncoupled. In this case Phya = but /3 ga i « Pb (since galaxies are 
predominantly dark matter). Also in table g, we display growth rate limits on a = 2/3^ 
when both galaxies and Lyman a fluctuations are uncoupled, P gSL \ = Phya = 0. This would 
be realized, for instance, in dark sector chameleon models where the fifth forces affects 
only cold dark matter and the effective coupling is density dependent. When Phya = 0, 
the constraints are very similar to those for Phya = Pb for < 3 but become stronger 
for larger z*. This is primarily due to the large observational value of fhya( z = 3) which, 
alone, slightly prefers Ph ya > 0. 



Table 1: The typical la constraints on the 
coupling a — 2/3q DM (z < z*) of the scalar 
field to CDM. We have considered two situ- 
ations. The first one corresponds to a uni- 
versal coupling of the scalar field to all the 
species (3 ga \ — /3cdm- In the second case, 
we assume that clustered objects such as 
galaxies have no coupling to the scalar field 
Pgai = 0. We have analysed the role of the 
coupling of the scalar field using galaxy sur- 
veys and Lyman a results. 



5. Concluding Remarks 

In this paper we considered a scalar-tensor theory with one scalar degree of freedom, whose 
coupling to matter is not universal. The couplings of cold dark matter, baryons, neutrinos, 
etc. to the scalar field is not the same, which will affect the growth of perturbations in the 
different matter species. 

To study large scale structure formation in this setup we considered two non-relativistic 
fluids: a dominant species B (playing the role of cold dark matter) and a subdominant 
species A (representing baryonic matter). We have found analytic formulae for the growth 
rate of the perturbations in both species if the couplings are constant. Furthermore we 
considered the case of perturbations crossing the Compton wavelength of the scalar field. 
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In doing so, we assumed that the (effective) coupling of the scalar field to the different 
matter species changed suddenly in the past. We were able to find approximate formulae 
for the growth rate. 

We then discussed the phenomenology of the 
theory and showed that there are new observa- 
tional signatures predicted by the theory such as 
an anomalous growth of structures influenced by a 
time dependent growth index 7 for both Cold Dark 
Matter and baryons. Most importantly, there is 
a single slip parameter S K relating the two New- 
ton potentials, but several ways of extracting this 
slip parameter from data. Our analysis provides a 
template for the growth of structures in the con- 
text of a well-motivated theoretical framework, i.e. 

, , n . , , , , , Table 2: The typical la constraints on 

scalar-tensor theories where the scalar couples to ? 

. the coupling a — 2I3q DM (z < z*) of the 
different matter species with constant couplings. . „ , , , T , ,, 

r scalar field to CDM. In both situations 

Different cosmological experiments (such those which displayed here the coupling to bary0 
probe the ISW, weak lensing or the distribution of 
galaxies) would probe the slip parameter and the 
growth of structures, and therefore the forces in 
the dark sector. It would be interesting to study 
how well future planned cosmological experiments 
can probe the type of effects discussed in this pa- 
per. 



2dFGRS, Lya and VVDS Data 


Phya = 


Phya = Pgal = 


z% 


a 




a 


0.5 


< 0.52 


0.5 


<57.4 


1 


<0.31 


1 


<6.7 


2 


<0.24 


2 


<1.6 


3 


<0.22 


3 


<0.87 


4 


<0.22 


4 


<0.65 
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yons 

and hence Lymana fluctuations is taken 
to vanish (fl^ya = 0). We have then con- 
sidered two situations. The first one cor- 
responds to a universal coupling of the 
scalar field to all dark matter /J ga i = 
/3cdm ■ In the second case, we assume that 
clustered objects such as galaxies have no 
coupling to the scalar field (i ga \ — 0. We 
have analysed the role of the coupling of 
the scalar field using galaxy surveys and 
Lyman a results. 
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